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Abstract. The nuclear symmetry energy at densities above saturation density (po ~ 
0.16/m~ 3 ) is poorly constrained theoretically and very few relevant experimental data 
exist. Its study is possible through Heavy Ion Collisions (HIC) at energies E/A > 200 
MeV, particularly with beams of neutron-rich radioactive nuclei. The energy range 
implies that the momentum dependence of the isospin fields, i.e. the difference of 
the effective masses on protons and neutrons, also has to be investigated before a 
safe constraint on E sym {p) is possible. We discuss the several observables which have 
been suggested, like n/p emission and their collective flows and the ratio of meson 
yields with different isospin projection, tt~/tt + and K°/K + . We point out several 
physical mechanisms that should be included in the theoretical models to allow a 
direct comparison to the more precise experiments which will be able to distinguish the 
isospin projection of the detected particles: CSR/Lanzhou, FAIR/GSI, RIBF/RIKEN, 
FRIB/MSU. 
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1. Introduction 

The study of the equation of state (EOS) of nuclear matter has been one of the main 
challenges of nuclear physics in the last 30 years. The main strategy for its determination 
has been the investigation of the pletheora of data from heavy ion experiments in a 
wide range of beam energies E/A ~ 0.02 — 2 GeV. In particular, through the analyisis 
of collective flows of nucleons and of strange meson production it has been possible to 
constrain the EOS in a wide range of densities up to p ~ 2p Pfl E]- This endeavor 
has been performed mainly with respect to the properties of symmetric nuclear matter 
with an equal number of protons and neutrons. However, it is known already from 
the Bethe-Weizsacker mass formula that the difference in the neutron-proton content 
gives a contribution to the energy called the symmetry energy, i.e. the energy arising 
from the breaking of the neutron-proton symmetry (apart form the Coulomb energy). 
The symmetry energy E syrn can be defined from the following expansion of the energy 
density e(p,p 3 ) = e(p) + pE sym (p)I 2 + 0(/ 4 ) + .., with p = p p + p n and p 3 = p p - p n 
being the total and the isospin densities and I = pzj p the isospin or asymmetry. 

The predictions on E sym (p) based on the existing many-body techniques diverge 
rather widely [2j [3] both at sub-saturation and supra-saturation densities, as shown in 
Fig.l (left) for several relativistic and non-relativistic models. Therefore, while E sym is 
fairly well constrained around saturation density from nuclear structure its behavior at 
lower and at higher density is poorly known [31 0]. 

The study of the E sym below po has been intensively pursued in the last ten years 
looking at a large variety of observables at different reaction mechanisms from giant 
dipole resonances, pre-equilibrium emission, multi- and neck-fragmentation, isospin 
diffusion, etc. [3j H] . The present contribution instead focuses on the open issues for the 
determination of E sym (p) in the high density region which is relevant for neutron star 
observables, like maximum masses, the mass-radius relation and cooling [6J, for hybrid 
star structure and the transition to a deconfined phase [7], and for the formation of 
black holes. In view of the planned experiments at GSI (ASY-EOS), and of experiments 
at the future facilities CSR/Lanzhou, FAIR/GSI, RIBF/RIKEN, FRIB/MSU[8J several 
observables which are sensitive to E sym to the high density region [3l HJ E] have been 
suggested, in particular the n/p and 3 He/t particle ratios and their collective flows and 
the isospin ratios of meson production 7r~ /tt + and K°/K + . We critically review these 
issues in this contribution. 

2. The Relativistic Structure of E sym 

The determination of E sym at high density p > 1.5po needs to exploit HIC's at 
beam energies of E/A > 200MeV in order that the maximum density is sufficiently 
high. For the theoretical approaches this shifts the energy regime of interest to a 
range where relativistic effects can start to become significant. Hence relativistic, 
covariant approaches should be preferred; nonetheless relativistic extensions of non- 
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relativistic approaches represent a viable way even if they may mix effects of relativity 
and many-body effects[10j. In this section we briefly recall the relativistic structure of 
E sym in effective models based on nucleons interacting through mesonic fields. At the 
level presented here the structure of E sym is shared by the different models based on 
effective Lagrangians of Quantum-Hadro-Dynamics (QHD), as the Relativistic Mean 
Field Theory (RMFT) [Til [12], the Density Dependent Relativistic Hartree approach 
(DDRH)P], the Effective Field Theory (EFT) density functional [14J as well as the 
microscopic DBHF[5j [15]. Therefore the following discussion allows to set a common 
language that will be useful in the discussion of meson production in Sect. [U 

In QHD the effect of nuclear interactions results in in-medium modifications of the 
scalar effective masses and the energy-momentum four-vector 

M; p = M + S CT ±S 5 , fc£ = *" + ££±E£, (1) 

where the different self-energies £j are labelled by the mesons representative of the 
specific spin-isospin quantum numbers of field. The upper and lower signs refer to 
neutrons and protons, respectively. The isovector self-energies can be written as 

E s (n,p) = -fspss , ^p(n,p) = fpjs , (2) 

where p s3 = p sn — p sp and j% = j£ — j£ are the scalar isospin density and the isospin 
current, respectively. In the most simple models of RMFT the coupling vertices fi are 
constant while in more sosphisticated models, as e.g. the DDHF, they are density and 
even momentum dependent [16J. 




Figure 1. Left and middle panel: Symmetry energy as a function of density as 
predicted by different models from Rcf.[2j. Right panel: Scalar and vector isovector 
self-energies in NL-RMFT; the solid line is to highlight the deviation from linearity. 



The symmetry energy arising from such a structure can be written as [12 
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where f p , fs are the couplings of the nucleons to the effective isovector scalar and vector 
field. However already from the Hartree-Fock theory as well as from the more complete 
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DBHF approach it is clear that each spin-isospin channel receives contributions from all 
the meson-like fields and the meson label is only indicative of the channel[17j. Eq.([3]) 
can be strictly derived only in RMFT, but for our purpose the key point is that E sym 
arises from a competition of two fields, one scalar and one vector, that at p ~ 2p Q are 
of the order of ~ 100 MeV, as is shown in Fig. 1 (right) for a standard non-linear 
(NL)-RMFT model [12]. This behaviour is shared by all relativistic approaches, with 
the difference that for example in DDRH the coupling f p ^ depends on density, similarly 
in DBHF. 

We stress the distinction between a scalar and a vector field because their balance 
in a system in equilibrium is generally dynamically broken for an evolving system. In 
Ref.[10j it has been shown that such an effect simulates a stiffer symmetry energy in the 
build-up of collective flows. In fact, the difference of the force acting on a neutron and 
a proton moving with momentum p can be written, after some approximations, as 

dp p _ dpn ^ 
dr dr 

where 7 is the Lorentz factor E/m. We can see that dynamically the vector field is 
enhanced and the scalar field suppressed resulting in an effective E* larger than E^"^ 
as can be seen by comparing Eqs.Q and (jlj). Hence, the relativistic structure implies a 
modified relation between E sym and observables in HIC's. Of course, one can simulate 
this relativistic structure with the help of a specific momentum dependence, but this 
entangles many-body effects with relativistic ones. A first analysis has shown that such 
effects are significant for E/A > 1 GeV [TO] . 

3. Isospin Momentum Dependence 

The many-body physics involved is indeed more complex than the one sketched above. 
In fact, the RMFT class of models are derived from a zero-range approximation that 
makes the meson propagator independent on the momentum. However, the explicit 
momentum dependence (MD) originating in the finite range of the interaction is known 
to be non- negligible in both relativistic and non- relativistic approaches. In fact, it has 
been shown that the self-energies of DBHF cannot be reproduced by the simple structure 
of RMFT, but an explicit MD is necessary. This has lead to the development of a RMFT 
with derivate coupling [16J that is able to reproduce the correct energy dependence of 
the optical potential. In Fig. 2 (left) the experimental optical potential is compared to 
the predictions of different theoretical approaches. We can see that at high momenta k 
also the DBHF is not able to reproduce U opt (k). 

The momentum dependence of the field is usually parametrized in terms of the 
(Landau) effective in-medium masses, defined as 



m q 



m 



^ m dU a 



p dp 



(5) 



where we use the index q = p, n to underline that, in general, they are different for 
proton and neutrons in asymmetric nuclear matter. The isospin independent part of 
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the optical model has been studied since many years and its behavior has been fairly well 
determined, see Fig. 2 (left, filled diamonds). Instead, the isospin dependence of the 
optical model, usually discussed in terms of the Lane potential U^ane = (U p — U n )/2I, 
shown in Fig. 2 (right), is poorly known, and this only in a restricted energy range, 
k < lfm^ 1 , where some data with large uncertainties are available. 




Figure 2. Nucleoli optical potential in nuclear matter at po. On the left side DBHF 
calculations for symmetric nuclear matter are compared to the phenomenological 
models NL3 and D3C and to the phenomenological p- A scattering analysis (diamonds) . 
The right panel compares the iso-vector optical potential for several relativistic and 
non-relativistic models. Figure taken from Ref.[2]. 

We will see that all the observables under investigation are in principle affected by 
the MD of the mean fields. This represents a further obstacle for the determination of 
E S ym especially at high p, because the importance of MD is expected to increase with 
energy. However, the possibility to have access experimentally to several observables as 
a function of momentum in a wide range offers the real possibility to constrain such a 
dependence. Generally, all RMFT approaches give m* < m*, while in non-relativistic 
approaches generally m* > m* with some except ions [18]. Non-relativistic BHF[20] 
would indicate the last choice as the correct one, which means that U L ane decreases with 
k. In RMFT one finds the opposite trend but it lacks the effect of the finite range of 
the interaction. In DBHF both relativistic fields and the finite range effect are included, 
but there exists an ambiguity from the method used to project on the different Lorentz 
amplitudes [2]. One could judge that the microscopic approaches favour the relation 
m* > m*, but in any case it is mandatory to show that this is the case in HIC by means 
of a comparison with the available and forthcoming data. 

Investigations of the MD of the isovector interaction have been performed using 
a non-relativistic transport code (Boltzmann-Nordheim-Vlasov, BNV) with a Skyrme- 
like mean field modified to account for a different MD for protons and neutrons [18j. 
At E/A = 400 AMeV for central collisions one can reach densities of about 1.8po for 
a time interval of the order of 15-20 fm/c. More interestingly, it is possible to find a 
kinematical region (y < 0.3) where the n/p and 3 H/ 3 He ratios of emitted particles as a 
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Figure 3. n/p ratio for 132 Sn + 124 Sn at 400 AMeV (b=lfm, y < 0.3) vs. the the 
transverse momentum normalized to the projectile one for two choices of E sym (p) and 

U L anc • 



function of the transverse momentum depend mainly on the neutron-proton splitting of 
the effective masses and are nearly independent of the stiffness of E sym . In Fig. 3 one 
sees that the n/p ratio exhibits a strong dependence for p t > p pro j = 2/m~ 1 depending 
on the sign of the splitting of the effective masses. For m* < m* the difference between 
the neutron and proton mean fields increases with momentum and the more energetic 
neutrons feel an increasing repulsion resulting in a larger n/p ratio. The opposite is 
true for the case m* < m*. The key point of our result [19] is that especially at high 
p t the effect is almost independent on the stiffness of the symmetry energy allowing to 
disentangle the dependence on the momentum from that on the density. Instead, for 
lower beam energies and in the low momentum region the pre-equilibrium emission is 
found to depend significantly on both the stiffness of E sym and the isospin MD. 

4. Isopin effects on Meson Production 

The study of the symmetry energy at high density implies to investigate HIC's at higher 
energy where a system at sufficiently high density is created. On the other hand, then 
the energy becomes large enough to excite the lower mass mesonic states, namely pions 
and kaons. Therefore symmetry energy effects are transferred to the production of 
mesons. In Ref.[2T] it was suggested that the 7i~/ri + ratio is sensitive to E sym (p). It 
was demonstrated that the main effect consists in the emission of nucleons during the 
evolution depending on E sym leading to a different n/p content of the residual system, 
which in turn influences the formation of the different charge states of the pion. 

Even if this mechanism seems quite straightforward one may doubt its real 
effectiveness. Indeed, the idea of using pions to determine the symmetric part of the EOS 
has suffered from the fact that pions strongly interact with nucleons and are produced 
during the whole evolution of the collision system making it difficult to associate their 
production to a specific density reached during the collision. In that context it was 
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Figure 4. Time evolution of the A ±,0 ' ++ resonances and pions n^' (left), and of 
kaons (K + '° (right) for a central (b = fm impact parameter) Au+Au collision at 1 
AGeV incident energy. Transport calculation using the NL, NLp, NLpS and DDF 
models for the iso-vector part of the nuclear EoS are shown. The inset shows the 
differential K°/K + ratio as a function of the kaon emission time. 

suggested by Aichelin and Ko[22] that kaons are a better probe of the EOS. The reason 
is manifold: Kaons have a higher threshold energy, hence they are produced only in the 
high density phase. Moreover, once produced they interact weakly with nucleons and 
their width with respect to the mass is quite small making a quasi-particle approximation 
more reliable. After nearly 20 years the effort to determine the symmetric EOS by kaon 
production has been successfull and is summarised in Ref . [23] . Following the same line 
of thinking the Catania group has suggested to investigate the K°/K + ratio as a better 
probe of the E sym at high density [2^1 [25] (the other isospin pair with anti-kaons K Q /K~ 
suffers from the strong coupling to the medium). 

Using a relativistic transport approach (Relativistic Boltzmann-Uhling-Uhlenbeck, 
RBUU[26]) we have analyzed pion and kaon production in central 197 Au+ 197 Au collisions 
in the 0.8 — 1.8 AGeV beam energy range, comparing models with the same "soft" 
EOS for symmetric matter and with different effective field choices for E sym . Fig. 4 
reports the temporal evolution of A ±,0,++ resonances, pions (tt ±, °) and kaons (K +, °) 
for central Au+Au collisions at lAGe^JJHJ [27] It is clear that, while the pion yield 
freezes out at times of the order of 50 fm/c, i.e. at the final stage of the reaction (and 
at low densities), kaon production occurs within the very early (compression) stage, 
and the yield saturates at around 15 fm/c, when the nucleon and A's densities reach 
their maximum value. From Fig. 4 we see that the pion multiplicities are moderately 
dependent on the isospin part of the nuclear mean field. However, a slight increase 
(decrease) in the tt~ (7r + ) multiplicity is observed when going from a NL-RMFT model 
without isovector fields (NL) to one with a p meson (NLp) and with p and 5 mesons 
(NLpd~). This trend is more pronounced for kaons, see the right panel, due to the high 
density selection of the source and the proximity to the production threshold. In Fig. 
5 (left) the pion and kaon rations from these calculations are shown as a function of 
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beam energy. The ratios decrease with beam energy, because the relative effects of 
mean fields and thresholds become less important. However, the greater sensitivitiy of 
the kaon ratios is seen clearly. 




Figure 5. Left: Excitation function of the tt~ /tt + and K°/K + ratios for Au + Au. 
RBUU results for different behavior of E«„, m (p)|25|. Right: Results for Ru + Ru in 
IBUU04 for a soft (x=l) and a stiff (x=0) -E sy m[2E] compared to the data from FOPI 
and a calculation with IQMD[29 . 

At present there are essentially no data for kaons, while there are some for pions 
and 7t~/tt + ratios from FOPI collaboration [29J. Within an isospin and momentum 
dependent transport model it has been shown that an agreement with data can be 
achieved only if a very soft E sym (p) is employed [28J. Such a finding appears in strong 
disagreement with other studies that exploit the elliptic flow to extract the slope of 
E sym around and above the saturation density [30]. However there are circumstantial 
reasons to be careful. Indeed, the physics involved in the in-medium particle production 
has many aspects and a comprehensive and self-consistent approach is necessary before 
extracting the isospin dependence of the interaction. 

The effect described in Ref. [28] essentially is due to the fact that a stiff E sym causes 
a neutron-rich emission of nucleons in the early stages of the reaction leaving the system 
too symmetric in isospin content to reproduce the 7r~/7r + ratio of FOPI. On the other 
hand, if one employs an E sym (p) that just above po decreases with density generating an 
isospin force that is attractive for the neutrons it is possible to get close to the data, see 
Fig. 5 (right). However, we note that it is mandatory to check if one can reproduce at 
the same time the n/p emission as a function of pt especially at the high p t relevant for 
pion production, since this represents the complementary signal for the soft E sym . This 
would be a first test of the claim for a very soft E sym and in case of failure it will provide 
evidence that there are other mechanism determining the in-medium meson production. 
Indeed, we already know that there are at least other three effects competing with 
the mean field effect on the n/p emission. These are the so-called "'threshold effect"' 
emphasized by the calculation of Ref. s [2HI25] , the momentum dependence of the isosp in- 
dependent cross sections, and the the isospin mass shift of pions due to the coupling to 
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A— hole excitation [21]. These will be discussed in the following. 
4-1. Isospin dependence of thresholds and spectral functions 

The "'threshold effect"' is due to the fact that the masses of nucleons and Delta's are 
modified in the medium. The unknown self-energies of the delta's are usually specified 
in terms of the neutron and proton ones by the use of Clebsch-Gordon coefficients for the 
isospin coupling of the A's to nucleons [2U [21] . These medium modifications are isospin 
dependent, Eq.(T5]). This influences the phase-space available for meson production in a 
nucleon-nucleon (NN) collision because it modifies the difference between the invariant 
energy in the entrance channel Sj„ with respect to the production threshold s t h- This 
effect is, of course, present in general for all meson productions, but for brevity we 
concentrate here on one inelastic channel: nn — > pA~ , mainly responsible for 7r~ 
production. From Eq.([T}) the invariant energy in the entrance channel and the threshold 
energy are given, respectively, by 



where m* N ,m\ are the isospin averaged values. The last equality for the threshold 
energy is valid due to the prescription for the Delta self-energies noted above, which 
leads to an exact compensation of the isospin-dependent parts, hence the threshold 
Sth is not modified by isospin dependent self-energies. In generally in a self-consistent 
many-body calculation higher order effects can destroy this exact balance. But this is 
not so important for our qualitative discussion, since there will always be a compensating 
effects in s t h- On the other hand, the energy available in the entrance channel, Sj n , is 
shifted in an explicitely isospin dependent way by the in-medium self-energy TP p +lls > 0. 
Especially, the vector self-energy gives a positive contribution to neutrons that increases 
the difference s in — s t h and hence increases the cross section of the inelastic process due 
to the opening up of the phase-space, expecially close to threshold. 

A similar modification but opposite in sign is present in s in — s t h for the pp — > n A ++ 
channel that therefore is suppressed by the isospin effect on the self-energies. Hence, 
due to the described threshold effect the ratio n~ /n + increases with the stiffness of E sym 
which is associated with a large S° + £5. This is at the origin of the result in Fig. 5 
(left). Of course, the RBUU calculation contains also the mean field effect on the isospin 
content but the final result appears to be dominated by the threshold effect. In fact, 
it is seen that tt~/tt + increases with the increasing stiffness of E sym and the effect is 
stronger closer to the threshold energy. 

The threshold effect described constitutes one physical effect that is present in the 
in-medium meson production and results in an opposite dependence on E sym relative 
to the mean field effect on the emission of neutrons and protons. On the other hand 
this does not mean that in Ref.[25] the effect is quantitatively described in the more 





(6) 
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realistic way. There are at least two other physical aspects that have to be considered. 
One is the fact that the self-energies in the implementation of RBUU is not explicitly 
momentum dependent as it is in DBHF and as the study of the optical potential 
shows, see Fig. 2 (right). Thus we return to the importance of the isospin momentum 
dependence discussed in Sect. El Because we neglect the dependence of the self-energies 
on momentum, the difference Si n — s t h in Eq.([6]) increases strongly with the momentum p. 
The problem is directly related to the strong energy dependence of the optical potential 
in RMFT in Fig. 2. Therefore it is likely that a more realistic calculation which includes 
a Y* Pt $(p,p) will reduce the effect seen in Fig. 5. In addition a fully consistent treatment 
should include an optical potential also for the pions. It can be envisaged that a stronger 
sensitivity to E(p,p) can be seen looking at n~ /n + as a function of the pion transverse 
energy. 

Another key point in the study of the isospin effects is the isospin dependence of the 
cross section for the inelastic processes. In fact the same medium effects that produce the 
threshold effect discussed above also modify the cross section of the pertinent inelastic 
processes. The problem of the self-consistency between the mean field and the cross 
section is known for many years, however, in the case under discussion it seems to be 
much more relevant. In fact one can simulate the effect of a repulsive mean field, e.g. 
on collective flow, by increasing the cross section. However, if one plays with the cross 
section for inelastic collisions it will manifest itself in the number of particles that are 
produced and measured. In this respect the calculations of pion and kaon yields of 
Ref. [25] take into account only one consequence of medium effects, namely the isopin 
dependent self-energies. It is likely that simultaneously taking into account the isospin 
dependence of the cross section can dampen the effect of self-energies on Si n — s t h- 

Another effect that is certainly present is the modification of the pion spectral 
function which in an asymmetric medium becomes isospin dependent. The effect has 
been pointed out in Ref. [31], see also the contribution of CM. Ko in these Proceedings. 
This work shows that the increased strength at low energies of the ir~ relative to the 
7r + spectral function in dense asymmetric nuclear matter modifies the expected thermal 
7t~/tt + ratio with respect to the one with vacuum masses. The effect is comparable 
in size to the mean field effect on n/p ratios of emitted nucleons, but again goes in 
the opposite direction: larger E sym implies larger ratios. Therefore it would 

be important to include such an effect in the transport models, even if it means to go 
beyond the simple quasi-particle approximation. From this point of view it would be 
desirable to have calculations with transport codes including off-shell cffccts[32j. 

Considering these various issues that are involved in the determination of the 
7r~//T + as well as the K°/K~ ratios the circumstantial evidence for a very soft E sym 
at high density should be considered only as a first step that opened the discussion. 
Instead, there is a need for theoretical improvements in the transport models to match 
the upcoming availability of experiments able to scrutiny more precisely the isopin 
dependence of observables in heavy ion collisions. 
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5. Outlook 



The determination of E sym (p) above the saturation density by means of HIC's is indeed 
at the beginning. With increasing energy the momentum dependence of the mean field 
determines the subsequent dynamics. Therefore the study of the E sym is entangled with 
the isospin effects on the momentum dependence. This affects not only the nucleon 
dynamics but also the isopin content of the meson production. We discussed several 
issues which from a theorethical point of view require further development of the present 
transport theories for the study of the isospin dependence of meson production. From 
the experimental point of view the possibility to look at observables especially as a 
function of momentum should allow to really constrain the E sym in the high density 
region. However a key point of the search for the high-density symmetry energy should 
be to look at several observables in the same reaction simultaneously which should help 
not to put under the carpet the important mechanisms at work in HIC dynamics. 
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